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Abstract: We report on further progress in understanding soft singularities of massless 
gauge theory scattering amplitudes. Recently, a set of equations was derived based on 
Sudakov factorization, constraining the soft anomalous dimension matrix of multi-leg scat- 
tering amplitudes to any loop order, and relating it to the cusp anomalous dimension. 
The minimal solution to these equations was shown to be a sum over color dipoles. Here 
we explore potential contributions to the soft anomalous dimension that go beyond the 
sum-over-dipoles formula. Such contributions are constrained by factorization and invari- 
ance under rescaling of parton momenta to be functions of conformally invariant cross 
ratios. Therefore, they must correlate the color and kinematic degrees of freedom of at 
least four hard partons, corresponding to gluon webs that connect four eikonal lines, which 
first appear at three loops. We analyze potential contributions, combining all available 
constraints, including Bose symmetry, the expected degree of transcendentality, and the 
singularity structure in the limit where two hard partons become collinear. We find that 
if the kinematic dependence is solely through products of logarithms of cross ratios, then 
at three loops there is a unique function that is consistent with all available constraints. 
If polylogarithms are allowed to appear as well, then at least two additional structures are 
consistent with the available constraints. 
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1. Introduction 

Understanding the structure of gauge theory scattering amplitudes is important from both 
the fundamental field-theoretic perspective, and the pragmatic one of collider phenomenol- 
ogy. Infrared singularities, in particular, open a window into the all-order structure of 
perturbation theory and the relation between the weak and strong coupling limits; at the 
same time, they provide the key to resummation of large logarithms in a variety of phe- 
nomenological applications. 

The study of infrared singularities in QCD amplitudes, which has a three-decade-long 
history [1-40], recently received a major boost [41-43]. The factorization properties of soft 
and collinear modes, also referred to as Sudakov factorization, were combined with the 
symmetry of soft-gluon interactions under rescaling of hard parton momenta, and were 
shown to constrain the structure of singularities of any massless gauge theory amplitude, 
to any loop order, and for a general number of colors N c . A remarkably simple structure 
emerges as the simplest solution to these constraints. All non-collinear soft singularities 
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are generated by an anomalous dimension matrix in color space [9, 10, 14, 15, 21, 28, 30, 
34,36]. In the simplest solution, this matrix takes the form of a sum over color dipoles, 
corresponding to pairwise interactions among hard partons. This interaction is governed 
by a single function of the strong coupling, the cusp anomalous dimension, 7x(«s)- The 
simplicity of this result is remarkable, especially given the complexity of multi-leg amplitude 
computations beyond tree level. The color dipole structure of soft singularities appears 
naturally at the one-loop order [24,28-30], where the interaction is genuinely of the form 
of a single gluon exchange between any two hard partons. The validity of this structure at 
two loops was not obvious a priori; it was discovered through the explicit computation of 
the anomalous dimension matrix [37, 38] . 

This remarkable simplicity is peculiar to the case of massless gauge theories: recent 
work [44-51] has shown that the two-loop matrix, when at least two colored legs are 
massive, is not proportional to the one-loop matrix, except in particular kinematic regions. 
In general, in the massive case, there are new contributions that correlate the color and 
momentum degrees of freedom of at least three partons, starting at two loops. These 
contributions vanish as 0(m 4 /s 2 ) in the small mass limit [49,50]. 

Given that all existing massless results are consistent with the sum-over-dipoles for- 
mula, it is tempting to conjecture that it gives the full answer [41-43,52,53]. As emphasized 
in Ref. [42], however, constraints based on Sudakov factorization and momentum rescaling 
alone are not sufficient to determine uniquely the form of the soft anomalous dimension. 
A logical possibility exists that further contributions will show up at the multi-loop level, 
which directly correlate the kinematic and color degrees of freedom of more than two hard 
partons. It is very interesting to establish whether these corrections exist, and, if they do 
not, to gain a complete understanding of the underlying reason. Beyond the significance 
of the soft singularities themselves, a complete understanding of their structure may shed 
light on the structure of the finite parts of scattering amplitudes. 

Ref. [42] showed that precisely two classes of contributions may appear as corrections 
to the sum-over-dipoles formula. The first class stems from the fact that the sum-over- 
dipoles formula provides a solution to the factorization-based constraints only if the cusp 
anomalous dimension, 7^(a s ), associated with a hard parton in representation i of the 
gauge group, obeys (a s ) = Cj7x(a s ), where 7ft- is universal and Cj is the quadratic 
Casimir of representation i, Ci = Ca or Cf for gluons and quarks, respectively. This 
property is referred to as 'Casimir scaling' henceforth. Casimir scaling holds through three 
loops [54,55]; an interesting open question is whether it holds at four loops and beyond [56]. 
At four loops, the quartic Casimir first appears in the color factors of diagrams for the cusp 
anomalous dimension. (In QCD, with gluons in the adjoint representation A, fermions in 
the fundamental representation F, and a Wilson line in representation R, the relevant 
quartic Casimirs are d a A bcd d a ^ cd and d F bcd d'^ cd , where d a ^ cd are totally symmetric tensors 
in the adjoint indices a,b,c,d.) However, Ref. [43] provided some arguments, based on 
factorization and collinear limits of multi-leg amplitudes, suggesting that Casimir scaling 
might actually hold at four loops. In the strong coupling limit, it is known to break down 
for M = 4 super- Yang- Mills theory in the large- N c limit [57] , at least when 7^ is computed 
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for Wilson lines in a special class of representations of the gauge group. 

The second class of corrections, the one on which we focus here, can occur even if 
the cusp anomalous dimension obeys Casimir scaling. In this case, the sum-over-dipoles 
formula solves a set of inhomogeneous linear differential equations, which follow from the 
constraints of Sudakov factorization and momentum rescalings. However, we can contem- 
plate adding solutions to the homogeneous differential equations, which are provided by 
arbitrary functions of conformally (and rescaling) invariant cross ratios built from the mo- 
menta of four hard partons [42]. Thus any additional terms must correlate directly the 
momenta, and colors, of four legs. Due to the non-Abelian exponentiation theorem [8,11,13] 
such contributions must originate in webs that connect four hard partons, which first ap- 
pear at three loops. From this perspective then, the absence of new correlations at two 
loops [37,38], or in three-loop diagrams involving matter fields [58], is not surprising, and 
it does not provide substantial new evidence in favor of the minimal, sum-over-dipoles so- 
lution. The first genuine test is from the matter-independent terms at three loops. At this 
order, purely gluonic webs may connect four hard partons, possibly inducing new types of 
soft singularities that correlate the color and kinematic variables of the four partons. 

The most recent step in addressing this issue was taken in Ref . [43] , in which an addi- 
tional strong constraint on the singularity structure of the amplitude was established, based 
on the properties of amplitudes as two partons become collinear. Recall that the primary 
object under consideration is the fixed-angle scattering amplitude, in which all ratios of 
kinematic invariants are taken to be of order unity. This fixed-angle limit is violated upon 
considering the special kinematic situation where two of the hard partons become collinear. 
An additional class of singularities, characterized by the vanishing invariant mass of the 
two partons, arises in this limit. The splitting amplitude is defined to capture this class 
of singularities. It relates an n-parton amplitude with two collinear partons to an (n — 1)- 
parton amplitude, where one of the legs carries the total momentum and color charge of 
the two collinear partons. The basic, universal property of the splitting amplitude is that 
it depends only on the momentum and color degrees of freedom of the collinear partons, 
and not on the rest of the process. 

Splitting amplitudes have been explicitly computed, or extracted from known scatter- 
ing amplitudes, at one [59-62] and two [63, 64] loops. A derivation of splitting- amplitude 
universality to all loop orders, based on unitarity, has been given in the large- N c limit [65]. 
The light-cone-gauge method for computing two-loop splitting amplitudes [63], in which 
only the two collinear legs and one off-shell parton appear, strongly suggests that the same 
all-orders universality extends to arbitrary color configurations, not just planar ones. 

Based on splitting-amplitude universality, Ref. [43] established additional constraints 
on the singularity structure of amplitudes. Using these constraints in conjunction with the 
Sudakov factorization constraints discussed above, that paper excluded any possible three- 
loop corrections depending linearly on logarithms of cross ratios. The final conclusion was, 
however, that more general functions of conformal cross ratios that vanish in all collinear 
limits could not be ruled out. 

In the present paper we re-examine the structure of soft singularities at three loops. We 
put together all available constraints, starting with the Sudakov factorization constraints 
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and Bose symmetry, and including the properties of the splitting amplitude and the ex- 
pected degree of transcendentality of the functions involved 1 . We make some plausible as- 
sumptions on the kinematic dependence, and consider all possible products of logarithms, 
and eventually also polylogarithms. We find that potential contributions beyond the sum- 
over-dipoles formula are still possible at three loops, but their functional form is severely 
constrained. 

The paper is organized as follows. We begin with three short sections in which we 
review the main relevant results of Refs. [42,43]. In Sec. 2 we briefly summarize the 
Sudakov factorization of the amplitude and the constraints imposed on the soft anomalous 
dimension matrix by rescaling invariance of Wilson lines. In Sec. 3 we present the sum- 
over-dipoles formula, the simplest possible solution to these constraints. In Sec. 4 we 
review the splitting amplitude constraint. The main part of our study is Sec. 5, in which 
we put together all available constraints and analyze the possible color and kinematic 
structures that may appear beyond the sum-over-dipoles formula. Most of the discussion 
is general, and applies to any loop order, but specific analysis is devoted to potential three- 
loop corrections. At the end of the section we make a few comments concerning four-loop 
corrections. Our discussion throughout Sec. 5 focuses on amplitudes involving four colored 
partons, plus any number of color-singlet particles. The generalization to the multi-parton 
case is presented in Sec. 6. Our conclusions are summarized in Sec. 7, while an appendix 
discusses the special case of four-parton scattering at three loops. 

2. Sudakov factorization and its consequences 

We summarize here the infrared and collinear factorization properties of fixed-angle scat- 
tering amplitudes A4 (pi/n, a s (/i 2 ), e) involving n massless partons, plus any number of 
color-singlet particles, evaluated in dimensional regularization with D = 4 — 2e. We refer 
the reader to Ref. [42] for technical details and operator definitions of the various functions 
involved. Multi-parton fixed-angle amplitudes can be expressed in terms of their color 
components Ml in a chosen basis in the vector space of available color structures for the 
scattering process at hand. All infrared and collinear singularities of Ml can be factor- 
ized [10,18,26,32,38,39,42] into jet functions Jj, one for each external leg i, multiplied by 
a (reduced) soft matrix Slm, 

M L {pi/^,a s {^ 2 ),e) = Slm (pij,a s (n 2 ),e) Hm { ~^ ' ~^^ , a s{fJ?),<^J 

x]W%^,a.(MV > ) - (2-1) 

1 Transcendentality here refers to the assignment of an additive integer r for each type of factor in a 
given term arising in an amplitude or an anomalous dimension: r = for rational functions, r = 1 for 
factors of it or single logarithms, hire; r = n for factors of ((n), In™ x or Li n (x), etc. [66]. We will provide 
more examples in Sec. 5.3. 
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leaving behind a vector of hard functions Hm, which are finite as e — > 0. A sum over 
M is implied. The hard momenta 2 pi defining the amplitude A4 are assumed to be light- 
like, pf = 0, while the rij are auxiliary vectors used to define the jets in a gauge- invariant 
way, and they are not light-like, nf / 0. The reduced soft matrix Slm can be computed 
from the expectation value of a product of eikonal lines, or Wilson lines, oriented along 
the hard parton momenta, dividing the result by n eikonal jet functions J, which remove 
collinear divergences and leave only singularities from soft, wide-angle virtual gluons. It is 
convenient to express the color structure of the soft matrix S in a basis-independent way, 
in terms of operators T", a = 1, 2, . . . , iV 2 — 1, representing the generators of S\J(N C ) acting 
on the color of parton i (i = 1, 2, . . . , n) [24]. 

The partonic (quark or gluon) jet function solves two evolution equations simul- 
tanously, one in the factorization scale /i and another in the kinematic variable (2pi-rii) 2 /nj 
(see e.g. Ref. [42]). The latter equation generalizes the evolution of the renormalization- 
group invariant form factor [22]. The resulting solution to these equations can be written 
as [53] 

Ji (^^"•O'V) = H Jt (l,a s ((W) >c ) exp | _ }_£ 2 ^ MaV)) 



T T /•(2 Pi .n i )Vn? dX 2 X / ^2 \ 1 

+ — J -» 4^^ (A ' e) ) 1 H(^^J + 2^ (as(A ' £)) 



(2.2) 



where Hj i is a finite coefficient function, and all singularities are generated by the exponent. 
The solution depends on just three anomalous dimensions, which are functions of the D- 
dimensional coupling alone: r yj i is the anomalous dimension of the quark or gluon field 
defining the jet (corresponding to the quantity 7* defined in Refs. [41,43]), while 7^ = 
2a s /ir + • • • and 5-g = a s /ir + ■ ■ ■ are, respectively, the cusp anomalous dimension and an 
additional eikonal anomalous dimension defined in Sec. 4.1 of Ref. [42]. In eq. (2.2) we 
have already assumed that the latter two quantities admit Casimir scaling, and we have 
factored out the quadratic Casimir operator Cj = Tj • TV 

Our main interest here is the reduced soft matrix S, which takes into account non- 
collinear soft radiation. It is defined entirely in terms of vacuum correlators of operators 
composed of semi- infinite Wilson lines (see e.g. Ref. [42]), and depends on the kinematic 
variables 

= {-fkjjj? __ [ft-/?/ e- 2i ^ 

Pn 2(/3 i -n i ) 2 2( /3j --n J -) 2 2 (ft • n t ) 2 2(ft • n,) 2 ' 1 " ) 

nf nj nf n 2 - 

which are invariant with respect to rescaling of all the Wilson line velocities fy. The fa are 
related to the external momenta by pf = {Q/y/2)P^ , where Q is a hard scale whose precise 
value will not be important here. The phases \j are defined by fa - (3j = —\fa- f3j\e~ mXij , 
where Ajj = 1 if i and j are both initial-state partons, or both final-state partons, and 



2 In our convention momentum conservation reads q + X^™=iP i = 0' wnere <7 is the recoil momentum 
carried by colorless particles. 
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Xij = otherwise. Note that the sign of the phase in e~ mXtj is determined by the +ie 
prescription for the Feynman propagator. 

The reduced soft matrix obeys the renormalization group equation 

H-^-S L m (Pij,a s ,e) = - ^Sln (pij,a s ,e) Y S NM {p ij ,a s ) . (2.4) 
a/x N 

The soft anomalous dimension matrix (pij,a s ), in turn, obeys the equation [42] 



S dl^p~- {Pij ' a ^ = J^K^s)S NM , V i, N, M , (2.5) 



found by considering a rescaling of the eikonal velocity The simplest solution of this 
equation is the sum-over-dipoles formula [42], 

rf ip a s ) = - - %(a s ) £ J] ln ftj T< • T, + - %(a s ) £ T, • T, . (2.6) 

i=l j't^i " i=l 

In this expression the dependence on the scale p appears exclusively through the argument 
of the D-dimensional coupling in and <5jj. Therefore eq. (2.4) is easily integrated to give 
the corresponding formula for the reduced soft matrix S, 



— f 1 d\ 2 \ 1 - n 

Sdip (Pij, a s , e) = exp j - - J — - <%(a s (A 2 , e)) ^ Tj • T, 

1 n 
--7x(a s (A 2 ,e)) ^^ln^T.-T 



(2.7) 



i=l j^i 



Equation (2.6) satisfies the constraints (2.5) if and only if the cusp anomalous dimension 
admits Casimir scaling, namely 7^(a s ) = Cijx( a s), with independent of the color 
representation of parton i. In this paper we shall assume that this is the case, postponing 
to future work the analysis of how higher-order Casimir contributions to jk would affect 
the soft anomalous dimension matrix (the starting point for such an analysis is eq. (5.5) 
of Ref. [42]). 

Even under the assumption of Casimir scaling for 7^, eq. (2.7) may not be the full 
result for S, because may receive additional corrections A 5 going beyond the sum-over- 
dipoles ansatz. In this case the full anomalous dimension can be written as a sum, 

T 3 ( Ptj ,a s ) = rf ip ( Pij ,a s ) + ( Ptj ,a s ) . (2.8) 

Here A 5 is a matrix in color space, which is constrained to satisfy the homogeneous differ- 
ential equation 

This equation is solved by any function of conformally invariant cross ratios of the form 

_ fk • Pj fh • Pi , oin . 
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which are related to the kinematic variables p^ in eq. (2.3), and to the momenta pi, by 

1/2 



Pijkl 



Pij Pkl 
Pik Pjl 



Pi ■ Pj Pk ■ Pi 
Pi ■ Pk Pj ■ Pi 



Pi ■ Pj Pk ■ Pi 



Pi ■ Pk Pj ■ Pi 



(2.11) 



Each leg that appears in p^ki does so once in the numerator and once in the denominator, 
thus cancelling in the combination of derivatives in eq. (2.9). Hence we define 



A 5 (pij,a s ) = A (pijki,a s ) . 



(2.12) 



Any additional correction A must introduce new correlations between at least four partons 
into the reduced soft function. Such additional corrections are known not to appear at two 
loops [37,38], as expected from the fact that two-loop webs can correlate at most three 
hard partons. By the same token, they cannot show up in matter-dependent diagrams at 
three loops, as verified explicitly in Ref. [58] . On the other hand, they might be generated 
at three loops by purely gluonic diagrams, such as the one shown in fig. 1. 




Figure 1: A purely gluonic diagram connecting the four hard partons labeled i,j,k,l, which may 
contribute to the soft anomalous dimension matrix at three loops. It correlates the colors of the 
four partons via the operator TfT^TgTf / ade / c&e . 



The main purpose of the present paper is to examine all available constraints on the 
soft anomalous dimension matrix, and check whether they are sufficient to rule out a non- 
vanishing A at three loops. We will show that, despite the powerful constraints available, 
corrections to the sum-over-dipoles formula may indeed appear at this order. In the case 
of purely logarithmic functions of the cross ratios Pijkl, we find a unique solution to all the 
constraints. Allowing also for the appearance of polylogarithms of a single variable, there 
are at least two additional solutions. 

3. Minimal ansatz for the singularities of the amplitude 

The factorization formula (2.1) has the attractive property that each of the singular factors 
is defined in a gauge- invariant way. It requires the introduction of the auxiliary vectors rii, 
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which have been very useful [42] in revealing the properties of the soft anomalous dimension. 
At the end of the day, however, the singularities of the amplitude M. cannot depend on 
these auxiliary vectors, but only on the kinematic invariants built out of external parton 
momenta. Indeed, as discussed below, the cancellation of the dependence of the singular 
terms on the vectors rti can be explicitly performed, and one can write the factorization of 
the amplitude in a more compact form: 

M^a s {i?)^ =z(ji-,a s (n}),e^J U (jj, ^l,a s (p 2 ),e^j , (3.1) 

as used in Refs. [41,43]. Here the (matrix-valued) Z factor absorbs all the infrared (soft and 
collinear) singularities, while the hard function % is finite as e — > 0. We distinguish between 
two scales, the renormalization scale which is present in the renormalized amplitude M. 
on the left-hand side of eq. (3.1), and [if, a factorization scale that is introduced through the 
Z factor. The function % (a vector in color space) plays the role of H in the factorization 
formula (2.1), but differs from it by being independent of the auxiliary vectors nj. 

Sudakov factorization implies that the Z matrix is renormalized multiplicatively. We 
can then define the anomalous dimension matrix T, corresponding to Z, by 

d ry ( Pi f 2 A \ ry ( Pi , 2\ 



(2L a .{fi),^ = -z(£L,a,(M}),e) T (^a.^. (3.2) 



Note that the matrix V is finite, but it can depend implicitly on e when evaluated as a 
function of the D-dimensional running coupling; it will then generate the infrared poles of 
Z, as usual, through integration over the scale. 

The sum-over-dipoles ansatz for V s , eq. (2.6), implies an analogous formula for T. 
In order to see it, one may use the factorization formula (2.1), substitute in eqs. (2.2) 
and (2.7), use color conservation, Ylj^a^j = ~ ^ii an d apply the identity 



(2pi ■ rii) 2 \ ( (2pj • n 3 



>2 



ln(---J + ln^-^] +ln 



/ \ 

(IA-&I e-^) 2 

\ nf n) j 



21n(2 \ Pi - Pj \ e- i7rAi 



(3.3) 

Note also that the poles associated with <%(a s ) cancel out between the soft and jet func- 
tions. In this way, one arrives at the sum-over-dipoles ansatz for T, 



r diP (f^(A 2 )) =-\ik K(A 2 )) EE ln 



=1 3& 



(3.4) 



+ E7J l («s(A 2 )) 



i=i 



The Z matrix which solves eq. (3.2) can be written in terms of the sum-over-dipoles 
ansatz (3.4) as an exponential, in a form similar to eq. (2.7). However, S^ip has only 
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simple poles in e in the exponent, while the integration of T^ip over the scale A of the 
D-dimensional running coupling will generate double (soft-collinear) poles within Z, inher- 
ited from the jet functions in eq. (2.2), because of the explicit dependence of Tdi p on the 
logarithm of the scale A. 

If a non-trivial correction A appears in the reduced soft function (2.8), then the full 
anomalous dimension is 

r (y,",(A 2 )) = T dip (p«,(A 2 )) + A { Pljkh a s (\ 2 )) . (3.5) 

In terms of this function the solution of eq. (3.2) takes the form 



z(^,a.^), £ ) = P ^p{-^" 2/ ^ r (p^(A 2 , e ))}, 0.6) 



where P stands for path-ordering: the order of the color matrices after expanding the 
exponential coincides with the ordering in the scale A. We emphasize that path ordering 
is only necessary in eq. (3.6) if A / and [A, Tdip] / 0. Indeed, the ansatz (3.4) has 
the property that the scale-dependence associated with non-trivial color operators appears 
through an overall factor, 7x(a s (A 2 )), so that color matrices T corresponding to different 
scales are proportional to each other, and obviously commute. This is no longer true for a 
generic A / 0, starting at a certain loop order /. In this case eq. (3.5) would generically 
be a sum of two non-commuting matrices, each of them having its own dependence on the 
coupling and thus on the scale A. Considering two scales Ai and A2, we would then have 
[r(Ai), r(A 2 )] ^ 0, and the order of the matrices in the expansion of eq. (3.6) would be 
dictated by the ordering of the scales. It should be noted, though, that the first loop order 
in Z that would be affected is order I + 1, because T starts at one loop, so that 

[r (AO , r (A 2 )] ~ [r« (Ai) , A« (A 2 )] = 0(a l s +1 ) . (3.7) 

The issue of ordering can thus be safely neglected at three loops, the first order at which 
a non- vanishing A can arise. 



4. The splitting-amplitude constraint 

Let us now consider the limit where two of the hard partons in the amplitude become 
collinear. Following Ref. [43], we shall see that this limit provides an additional constraint 
on the structure of A. The way we use this constraint in the next section will go beyond 
what was done in Ref. [43]; we will find explicit solutions satisfying the constraint (as well 
as other consistency conditions discussed in the next section). 

The Sudakov factorization described by eq. (2.1), and subsequently the singularity 
structure encoded in Z in eq. (3.1), apply to scattering amplitudes at fixed angles. All 
the invariants pi ■ pj are taken to be of the same order, much larger than the confinement 
scale A 2 . The limit in which two of the hard partons are taken collinear, e.g. p\ ■ p2 — > 0, 
is a singular limit, which we are now about to explore. In this limit, p\ — > zP and 
P2 — > (1 — z)P, where the longitudinal momentum fraction z obeys < z < 1 (for time-like 
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splitting). We will see, following Ref. [43], that there is a relation between the singularities 
that are associated with the splitting — the replacement of one parton by two collinear 
partons — and the singularities encoded in Z in eq. (3.1). 

It is useful for our derivation to make a clear distinction between the two scales [if 
and [i introduced in eq. (3.1). Let us first define the splitting amplitude, which relates the 
dimensionally-regularized amplitude for the scattering of n — 1 partons to the one for n 
partons, two of which are taken collinear. We may write 

X || 2 

M n (pi,p2,Pj;n,e) — > Sp(pi,p 2 ;n,e) M n -i (P,Pj-,fJ,,e) . (4.1) 

Here the two hard partons that become collinear are denoted by p\ and p 2 , and all the 
other momenta by pj, with j = 3,4, ... ,n. We have slightly modified our notation for 
simplicity: the number of colored partons involved in the scattering is indicated explicitly; 
the dependence of each factor on the running coupling is understood; finally, the matrix 
elements have dimensionful arguments (while in fact they depend on dimensionless ratios, 
as indicated in the previous sections). The splitting described by eq. (4.1) preserves the 
total momentum p\ +p 2 = P and the total color charge Ti + T2 = T. We assume eq. (4.1) 
to be valid in the collinear limit, up to corrections that must be finite as P 2 = 1p\ -p2 — > 0. 

The splitting matrix Sp encodes all singular contribution to the amplitude M n arising 
from the limit P 2 — > 0, and, crucially, it must depend only on the quantum numbers of 
the splitting partons. The matrix element M n -i, in contrast, is evaluated at P 2 = 0, and 
therefore it obeys Sudakov factorization, eq. (3.1), as applied to an (n— l)-parton amplitude. 
The operator Sp is designed to relate color matrices defined in the ra-parton color space 
to those defined in the (n — l)-parton space: it multiplies on its left the former and on 
its right the latter. Thus, the initial definition of Sp is not diagonal. Upon substituting 
T = Ti + T2, however, one can use the n-parton color space only. In this space Sp is 
diagonal; all of its dependence on Ti and T2 can be expressed in terms of the quadratic 
Casimirs, using 2 T 1 ■ T 2 = T 2 - T\ - T 2 ,. 

Because the fixed-angle factorization theorem in eq. (2.1) breaks down in the collinear 
limit, p\-p2 — > 0, we expect that some of the singularities captured by the splitting matrix 
Sp will arise from the hard functions %. Specifically, if the Z factor in eq. (3.1) is defined 
in a minimal scheme, % will contain all terms in Ai n with logarithmic singularities in p\ -p2 
associated with non-negative powers of e. We then define Sp n , in analogy with eq. (4.1), 
by the collinear behavior of the hard functions, 

x 11 2 

'Hn(pi,P2,Pj;v,Hf,e) — > Sp n (p 1 ,p 2 ;n,fif,e)Hn-i {P,pj; n, Hf,e) , (4.2) 

where all factors are finite as e — > 0. As was the case for eq. (4.1), eq. (4.2) is valid up 
to corrections that remain finite in the limit P 2 — > 0. Singularities in that limit are all 
contained in the splitting matrix Sp n , while the function H n -i is evaluated at P 2 = 0. 

Next, recall the definition of the Z factors in eq. (3.1) for both the n- and (n— l)-parton 
amplitudes. In the present notation, they read 

M n {pi,P2,Pj-,IJ',e) = Z n (p 1 ,p 2 ,p j ;nf,e) H n (pi,P2,Pj;v,IJ'f,e) , (4.3) 
M n -i (P,Pj;fJ,,e) = Z n _ x (P,pj;n f ,e) U n -\ (P,Pf, V, Vf, e) ■ (4.4) 
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Substituting eq. (4.4) into eq. (4.1) yields 



1||2 

M n (pi,p2,Pj;^,e) — > Sp{p 1 ,p 2 ;fi,e)Z n - 1 (P,p j ;fif,e) U n -\ (P,Pf, M, P>f, e) • (4.5) 
On the other hand, substituting eq. (4.2) into eq. (4.3) we get 

X II 2 

M n (pi,P2,Pj-,tJ-,€) — > Z n (p 1 ,p 2 ,p j ;nf,e) Sp n (p 1 ,p 2 ; fx, ///, e) H n -i (P,Pj-, fx, /x/, e) . 

(4.6) 

Comparing these two equations we immediately deduce the relation between the full split- 
ting matrix Sp, which is infrared divergent, and its infrared-finite counterpart Sp%, 

Spft(pi,f>2;/x,M/,e) = Z' 1 (p 1 ,p 2 ,p j ;nf,e) Sp(p 1 ,p 2 ; fx, e) Z n _i (P,p.,-; fx/, e) , (4.7) 

where Z n is understood to be evaluated in the collinear limit. This equation (cf. eq. (55) in 
Ref. [43]) is a non-trivial constraint on both Z and the splitting amplitude Sp, given that 
the left-hand side must be finite as e — > 0, and that the splitting amplitude depends only 
on the momenta and color variables of the splitting partons — not on other hard partons 
involved in the scattering process. 

To formulate these constraints, we take a logarithmic derivative of eq. (4.7), using the 
definition of T n and T n _i according to eq. (3.2). Using the fact that Sp(pi,p 2 ; fx, e) does 
not depend on /xj, it is straightforward to show that 

■Sp H {pi,P2;fx,tXf,e) = T n (pi,P2,Pj-,fJ>f) Sp n (p 1 ,p 2 ;n,Hf,e) 



- Sp w (pi,p 2 ;^,M/,e) r n-i (P,Pj-,fJ,f) , 

where, as above, (n — l)-parton matrices are evaluated in collinear kinematics (P 2 = 0), 
and corrections are finite in the collinear limit. Note that all the functions entering (4.8) 
are finite for e — > 0. Note also that we have adapted the T matrices to our current notation 
with dimensionful arguments; as before, the matrices involved acquire implicit e dependence 
when evaluated as functions of the D-dimensional coupling. 

Upon using the identification T = Ti + T2, the matrix T n _i can be promoted to 
operate on the n-parton color space. Once one does this, one immediately recognizes that 
the splitting matrix Sp commutes with the T matrices, as an immediate consequence of 
the fact that it can only depend on the color degrees of freedom of the partons involved in 
the splitting, i.e. Ti, T2 and T = Ti + T2, and it is therefore color diagonal. Therefore, 
we can rewrite eq. (4.10) as an evolution equation for the splitting amplitude: 

-rr — s P«(pi>f2;^/,e) = r Sp (pi,P2;/v) s P«(pi>P2;^/,e) , (4.9) 

where we defined 

rsp(pi,f>2;/v) = r n {pi,P2,Pj\vf) - r n _i (p,pj;n f ) . (4.10) 
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We may now solve eq. (4.9) for the fif dependence of Sp n , with the result 



Sp w (pi,P2;/^/; e ) = Sp^\p 1 ,p 2 ;n,e) exp 



1 /■"/ d\ 2 



The initial condition for evolution 



Sp-^Pi^P,^/ = /i,e) 



(4.11) 



(4.12) 



will, in general, still be singular as pi -p 2 - > 0, although it is finite as e — ?■ 0. We may, in any 
case, use eq. (4.11) by matching the ^-dependence in eq. (4.7), which yields an expression 
for the full splitting function Sp. We find 



Sp(pi,2?2;/i,e) = Sp^ (p 1 ,p 2 ;fi,e) exp 



1 r 2dx \ ( n 

2 7 "p-^P^ 1 '^ 25 ^ 



(4.13) 



While collinear singularities accompanied by non-negative powers of e are still present in 
the initial condition, all poles in e in the full splitting matrix arise from the integration 
over the scale of the D-dimensional running coupling in the exponent of eq. (4.13). 

The restricted kinematic dependence of Ts p , which generates the poles in the splitting 
function Sp, is sufficient to provide nontrivial constraints on the matrix A, as we will now 
see. Indeed, substituting eq. (3.5) into eq. (4.10) we obtain 



rs P (pi,P2;A) = rs Pj dip(pi,P2;A) + A„(py W ;A) - A n _i(py W ;A) 



(4.14) 



where 



rs P ,dip(pi,P2;A) =- -%(a s (A 2 )) 



In 



2 \pi -P2\ e 



-\tt\t_2 



A 2 



Ti-T 2 



- Ti • (Ti + T 2 ) \nz - T 2 ■ (T x + T 2 ) ln(l - z) 

+ 7Jl {a s {\ 2 )) + 7J2 (a s (A 2 )) - 7 Jp («s(A 2 )) . 



(4.15) 



Equation (4.15) is the result of substituting the sum-over-dipoles ansatz (3.4) for T n and 
r n _i. The terms in eq. (4.14) going beyond eq. (4.15) depend on conformally invariant cross 
ratios in the ra-parton and (n — l)-parton amplitudes, respectively. Their difference should 
conspire to depend only on the kinematic variables p\ and p 2 and on the color variables 
Ti and T 2 . In this way eq. (4.14) provides a non-trivial constraint on the structure of A, 
which we will implement in Sec. 5.4. 



5. Constraining corrections to the sum-over-dipoles formula 
5.1 Functions of conformally invariant cross ratios 

Our task here is to analyze potential contributions of the form A (Pijkh&s) to the soft 
singularities of any n-leg amplitude. Our starting point is the fact that these contribu- 
tions must be written as functions of conformally invariant cross ratios of the form (2.10). 
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Because we are dealing with renormalizable theories in four dimensions, we do not expect 
Z to contain power-law dependence on the kinematic variables; instead the dependence 
should be "slow", i.e. logarithmic in the arguments, through variables of the form 



Eventually, at high enough order, dependence on Pijki through polylogarithms and har- 
monic polylogarithms might arise. We will not assume here that A is linear in the variables 
Lijki- We will allow logarithms of different cross ratios to appear in a product, raised to 
various powers, and this will be a key to finding solutions consistent with the collinear 
limits. Subsequently, we will examine how further solutions may arise if polylogarithmic 
dependence is allowed. 

A further motivation to consider a general logarithmic dependence through the vari- 
ables in eq. (5.1) is provided by the collinear limits, which can take certain cross ratios Pijki 
to 0, 1, or oo, corresponding to physical limits where logarithmic divergences in A will be 
possible. Other values of the cross ratios, on the other hand, should not cause (unphysical) 
singularities in A(py^). This fact limits the acceptable functional forms. For example, 
in specifying a logarithmic functional dependence to be through eq. (5.1), we explicitly 
exclude the form ln(c + Piju) for general 3 constant c. Such a shift in the argument of 
the logarithm would generate unphysical singularities at Pijki = — c, and would also lead 
to complicated symmetry properties under parton exchange, which would make it diffi- 
cult to accomodate Bose symmetry. We will thus focus our initial analysis on kinematic 
dependence through the variables L^ki- Although it seems less natural, polylogarithmic 
dependence on p^i cannot be altogether ruled out, and will be considered in the context 
of the three-loop analysis in Sec. 5.5. 

The fact that the variables (5.1) involve the momenta of four partons, points to their 
origin in webs that connect (at least) four of the hard partons in the process, exemplified 
by fig. 1. The appearance of such terms in the exponent, as a correction to the sum- 
over-dipoles formula, implies, through the non-Abelian exponentiation theorem, that they 
cannot be reduced to sums of independent webs connecting just two or three partons, nei- 
ther diagrammatically nor algebraically. Indeed, for amplitudes composed of just three 
partons the sum-over-dipoles formula is exact [42]. Similarly, because two- loop webs can 
connect at most three different partons, conformally invariant cross ratios cannot be formed. 
Consequently, at two loops there are no corrections to the sum-over-dipoles formula, inde- 
pendently of the number of legs. Thus, the first non-trivial corrections can appear at three 
loops, and if they appear, they are directly related to webs that connect four partons. For 
the remainder of this section, therefore, we will focus on corrections to the sum-over-dipoles 
formula that arise from webs connecting precisely four partons, although other partons or 
colorless particles can be present in the full amplitude. Our conclusions are fully general 
at three loops, as discussed in Sec. 6, because at that order no web can connect more than 
four partons. 

3 In Sec. 5.6 we will briefly consider the possibility of including a dependence of the form ln(l — p). 




(5.1) 
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We begin by observing that, independently of the loop order at which four-parton 
corrections appear, their color factor must involve at least one color generator corresponding 
to each of the four partons involved. For example, the simplest structure a term in A can 
have in color space is 

MPijkl) = h abcd T" T$ T c k Tf &f n ( Pm ) , (5.2) 

where h abcd is some color tensor built out of structure constants corresponding to the 
internal vertices in the web that connects the four partons (i,j,k,l) to each other. Note 
that h abcd may receive contributions from several different webs at a given order, and 
furthermore, for a given h abcd , the kinematic coefficient A^ m (pijki) can receive corrections 
from higher-order webs. In what follows, we will not display the dependence on the coupling 
of the kinematic factors, because it does not affect our arguments. As we will see, symmetry 
arguments will, in general, force us to consider sums of terms of the form (5.2), with different 
color tensors h abcd associated with different kinematic factors. 

More generally, at sufficiently high orders, there can be other types of contributions in 
which each Wilson line in the soft function is attached to more than one gluon, and hence 
to more than one index in a color tensor. Such corrections will be sums of terms of the 
form 



^(Pijkl) =Af n (p ijk l) /i i'-.am 1) 6i,...,6 m2) ci 1 ... 1 c msi di 1 ...,d m4 

\ i i ^33 j k k k I I I 



(5-3) 



where ()+ indicates symmetrization with respect to all the indices corresponding to a 
given parton. Note that generators carrying indices of different partons commute, while 
the antisymmetric components have been excluded from eq. (5.3), because they reduce, 
via the commutation relation [Tf , T b ] = if abc T^, to shorter strings 4 . In the following 
subsections, we will focus on (combinations of) color structures of the form (5.2), and we 
will not consider further the more general case of eq. (5.3), which, in any case, can only 
arise at or beyond four loops. 

5.2 Bose symmetry 

The Wilson lines defining the reduced soft matrix are effectively scalars, as the spin- 
dependent parts have been stripped off and absorbed in the jet functions. Consequently, 
the matrices V and A should admit Bose symmetry and be invariant under the exchange of 
any pair of hard partons. Because A depends on color and kinematic variables, this sym- 
metry implies correlation between color and kinematics. In particular, considering a term 
of the form (5.2), the symmetry properties of h abcd under permutations of the indices a, b, 
c and d must be mirrored in the symmetry properties of the kinematic factor A\ m (pijj c i) 
under permutations of the corresponding momenta pi, pj, pk and p\. The requirement of 



4 One can make a stronger statement for Wilson lines in the fundamental representation of the gauge 
group. In that case the symmetric combination in eq. (5.3) can also be further reduced, using the identity 
{t a ,t(,} = jrS a b + d a bctc, so that the generic correction in eq. (5.3) turns into a combination of terms of 
the form (5.2). We are not aware of generalizations of this possibility to arbitrary representations. 
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Bose symmetry will lead us to express A as a sum of terms, each having color and kinematic 
factors with a definite symmetry under some (or all) permutations. 

Because we are considering corrections arising from four-parton webs, we need to 
analyze the symmetry properties under particle exchanges of the ratios pijki that can be 
constructed with four partons. There are 24 different cross ratios of this type, corresponding 
to the number of elements of the permutation group acting on four objects, S4. However, 
a Z2 x Z2 subgroup of S4 leaves each p^ki (and hence each L^ki) invariant. Indeed, one 
readily verifies that 

Pijkl = Pjilk = Pklij = Plkji ■ (5.4) 

The subgroup Z2 X Z2 is an invariant subgroup of S4. Thus, we may use it to fix one of 
the indices, say i, in p^ki- This leaves six cross ratios, transforming under the permutation 
group of three objects, S3 ~ S±/{Z2 X Z2). 



P 1234 




P 1324 1 < >2 1< > 3 



Figure 2: Symmetry properties of conformally invariant cross ratios. Each of the two triangles in 
the left-hand figure connects three cross ratios. The cross ratios associated with the two triangles are 
related by inversion, and one moves between the two triangles with odd permutations of momentum 
labels. The right-hand figure shows the resulting antisymmetry of the logarithms of the three 
conformally invariant cross ratios under permutations. The three variables L1234, £1423 and L1342 
transform into one another under permutations, up to an overall minus sign. For example, under 
the permutation 1 H 2, we have L1234 — > —£1423, £1423 -> —£1234, and £1342 — > —£1342- 



The permutation properties of the remaining six cross ratios are displayed graphically 
in fig. 2, where we made the identifications {i,j,k,l} — > {1,2,3,4} for simplicity. The 
analysis can be further simplified by noting that odd permutations in S4 merely invert 
Pijki, so that, for example, 

Pijkl = — * Lijki = — Likji ■ (5-5) 

Pikjl 

This inversion corresponds to moving across to the diametrically opposite point in the 
left-hand plot in fig. 2. We conclude that there are only three different cross ratios (corre- 
sponding to the cyclic permutations of {j, k,l} associated with S3/Z2 — Z3), namely Pijki, 
pujk and pikij- They correspond to triangles in fig. 2. Finally, the logarithms of the three 
cross ratios are linearly dependent, summing to zero: 

Lijki + Lnjk + Likij = . (5.6) 
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These symmetry properties lead us to consider for A4 111 in eq. (5.2) the general form 
Af n ( Pljkl ) = (L 12 34) hl (Lu 23 ) h2 (L 1U2 ) h3 , (5.7) 

where we have adopted the labeling of hard partons by {1, 2, 3, 4} as in fig. 2. Here the hi 
are non-negative integers, and eq. (5.6) has not yet been taken into account. Our general 
strategy will be to construct linear combinations of the monomials in eq. (5.7) designed to 
match the symmetries of the available color tensors, h ahcd in eq. (5.2). Such combinations 
can be constructed for general h{. As we shall see, however, transcendentality constraints 
restrict the integers hi to be small at low loop orders. In the three-loop case, this will 
suffice to eliminate all solutions to the constraints, except for a single function. 

We begin by noting that the antisymmetry of L1234 under the permutation lf*4 (or 
under 2 -H- 3, see fig. 2) is mirrored by the antisymmetry of the color factor h abcd T^T^TgT^ 
if the tensor h abcd = f ade f cbe ^ where f ade are the usual, fully antisymmetric SU(iV c ) struc- 
ture constants. The same is obviously true for any odd power of £1234, while in the case of 
even powers, an appropriate type of color tensor is h abcd = d ade d cbe , where d ade are the fully 
symmetric SU(iV c ) tensors. Fig. 2, however, shows that under other permutations, the dif- 
ferent cross ratios transform into one another. Therefore, if we are to write a function with 
a definite symmetry under all permutations, it must be a function of all three variables. 
Specifically, in order for a term of the form (5.7) to have, by itself, a definite symmetry 
under all permutations, the powers h±, h 2 and /13, must all be equal. Alternatively, one 
can consider a linear combination of several terms of the form (5.7), yielding together a 
function of the kinematic variables with definite symmetry. In this respect it is useful to 
keep in mind that the sum of the three logarithms (all with a single power) is identically 
zero, by eq. (5.6). 

Let us now construct the different structures that realize Bose symmetry, by considering 
linear combinations of terms of the form of eq. (5.2), with A4 111 given by eq. (5.7). We 
consider first three examples, where the logarithms are raised to a single power h. As 
we will see, none of these examples will satisfy all the constraints; they are useful however 
for illustrating the available structures. 

a) We first consider simply setting h\ = h 2 = /13 in eq. (5.7), obtaining 

^(Pijkl) = h abcd TJT2T3T4 L1234 -^1423 ^1342 • (5.8) 



For odd h the color tensor h abcd must be completely antisymmetric in the four indices, 
while for even h it must be completely symmetric. We anticipate that odd h is ruled 
out, because completely antisymmetric four- index invariant tensors do not exist for 
simple Lie groups [67]. Furthermore, while symmetric tensors do exist, eq. (5.8) is 
ruled out at three loops, because from fig. 1 it is clear that only h abcd = f ade f che ( Q r 
permutations thereof) can arise in Feynman diagrams at this order. 

b) Our second example is 

A 4 ( Pijkl ) = T?T|T§t2 \f ade f cbe L h 12U + f cae f dbe L h U23 + f bae f cde L h lU2 \ , (5.9) 
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where h must be odd. Alternatively, each f ade may be replaced by the fully symmetric 
SU(iV c ) tensor d ade , and then h must be even. In eq. (5.9) each term has a definite 
symmetry only with respect to certain permutations, but the three terms transform 
into one another in such a way that their sum admits full Bose symmetry. We will 
see shortly that the structure in eq. (5.9) does not satisfy the collinear constraints. 

c) Finally, one may consider the case where two of the three logarithms appear together 
in a product, raised to some power h, 

A l „ \ rriar-pbrpcrpd jabe icde i j j \h 

AKPijkl) = J-l - L 2- L 3 1 4 « « \ L 12U -^1423 ) 

n (5-10) 

+ d dae d cbe (L U23 L 1U2 ) h + d cae d bde (L 1U2 L 12U ) h . 

Once again, we observe that these color tensors cannot arise in three-loop webs. 
Furthermore, as we will see, eq. (5.10), at any loop order, fails to satisfy the collinear 
constraints. 

We are led to consider more general structures, using eqs. (5.2) and (5.7) with arbitrary 
integers h{. As announced, we will satisfy Bose symmetry by constructing polynomial 
kinematical factors mimicking the symmetry of the available color tensors. One may write 
for example 



MPijkl) = T?T£T§T2 

x 



rade fcbe j ti\ I r/12 r n 3 ( _-> \/ii+ft 2 +l3 T " 2 T 

J J -^1234 ^1423 -^1342 V 1 ) ^1342 ^14 

1 pcaefdbe T h\ ( r/12 7-/13 I -I\hi+h2+h3 r^2 t 

+ / / ^1423 ^1342-^1234 ~~ ^1234 ^ 

1 fbaefcde t hi ( j hi 7-/13 / -1 \hi+h 2 +h 3 jh 2 j 

-r J J -^1342 ^1234 "^1423 V L ) ^1423 ^ 



hz 

1234 -^1342 



(5.11) 



h 3 

1423 "^1234 



where hi, h 2 and /13 can be any non-negative integers. The first line is invariant, for 
example, under the permutation 1 o 4 (when applied to both kinematics and color), the 
second line is invariant under 1 <H> 3, and the third is invariant under 1 -H- 2. The other 
exchange symmetries are realized by the transformation of two lines into one another. 
For example, under 1 4 the second line transforms into the third and vice versa. In 
eq. (5.11) the color and kinematic factors in each line are separately antisymmetric under 
the corresponding permutation. Note that eq. (5.9) corresponds to the special case where 
hi in eq. (5.11) is odd, while h 2 = /13 = 0. 

One can also construct an alternative Bose symmetrization using the symmetric com- 
bination, 



x 



■ade icbe jh\ I jh 2 7-/13 . ( _-i\hi+h 2 +h 3 jh 2 7-/13 

u -^1234 \^1423 "^1342 "•" V 1 ) -^1342 "^1423 

1 jcae jdbe jh\ ( j h 2 j h 3 . / 1 \hi+h2+h 3 j h 2 7-/13 \ 

a -^ 142 3 ^1342-^1234 V L ) -^1234 ^1342^ 

1 jbaeicde j hi ( 7-/12 7-/13 1 I 1 \hi+h 2 +h 3 jh 2 7-/13 \ 

-I- a a -^1342 ^1234-^1423 ~r I L ) -^1423 ^1234^ 



(5.12) 
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Note that eq. (5.10) is reproduced by setting hi = and hi = /13 = h in eq. (5.12). 

Eqs. (5.11) and (5.12) both yield non-trivial functions for both even and odd powers hi, 
with the following exceptions: For even hi, eq. (5.11) becomes identically zero if hi = /13; 
similarly, for odd hi eq. (5.12) becomes identically zero if hi = h%. It is interesting to note 
that eq. (5.8) with odd h cannot be obtained as a special case of eq. (5.11). Indeed, by 
choosing hi = hi = /13 one obtains the correct kinematic dependence, but then the color 
structure factors out and vanishes by the Jacobi identity, 

j^abcd jade jcbe _j_ jcae jdbe _j_ jbae jcde q (5 13) 

In contrast, for even h, eq. (5.8) can be obtained as a special case of eq. (5.12), setting 

h abcd = d ade d cbe + ^cae^dbe + jbaejcde ^ ^ ^ 

which is totally symmetric, as required. This is expected from the general properties of 
symmetric and antisymmetric invariant tensors for simple Lie algebras [67]. 

At any fixed number of loops /, the total power of the logarithms in eqs. (5.11) 
and (5.12), /i to t = hi + hi + /13, will play an important role. Indeed, h tot is the degree 
of transcendent ality of the function A4, as defined in the Introduction, and it is bounded 
from above by the maximal allowed transcendentality of the anomalous dimension matrix 
at I loops, as described in Sec. 5.3. We expect then that at / loops there will be a finite 
number of sets of integers hi satisfying the available constraints. The most general solution 
for the correction term A4 will then be given by a linear combination of symmetric and 
antisymmetric polynomials such as those given in eqs. (5.11) and (5.12), with all allowed 
choices of hi. 

Such combinations include also contributions related to higher-order Casimir opera- 
tors. Indeed, summing over permutations of {hi,hi,h^} in the symmetric version of A4, 
eq. (5.12), one finds a completely symmetric kinematic factor, multiplying a color tensor 
which is directly related to the quartic Casimir operator (with a suitable choice of basis in 
the space of symmetric tensors over the Lie algebra [67]), 

A 4 (Pijfci) = T?T|T§Tf x \d ade d cbe + d cae d dbe + d hae d cde 

jh\ jh.2 r/13 1 rhi jhi jhz , j-hi r/12 r/13 {r i r:\ 

-^1234 "^1423 "^1342 "^1423 -^1342 "^1234 t -^1342 "^1234 "^1423 [0.10) 

I ( 1 \hi+h2+h3 ( t hi j\i2 t h 3 , T h\ jh 2 t h s , jh\ jh 2 t h 3 \ 
I L ) ^1234 "^1342 -^1423 "^1423 "^1234 "^1342 -^1342 -^1423 -^1234 J ■ 

For even htot = hi + hi + h^ this function is always non-trivial, while for odd h tot it is only 
non-trivial if all three powers hi are different. We note once again that, due to the Jacobi 
identity (5.13), eq. (5.15) does not have an analog involving the antisymmetric structure 
constants. 

5.3 Maximal transcendentality 

Our next observation is that, at a given loop order, the total power of the logarithms, /itot> 
cannot be arbitrarily high. It is well known (although not proven mathematically) that the 
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maximal transcendentality r max of the coefficient of the 1/e pole in an Z-loop amplitude 
(including k = 0) is T max — 21 — fc. If a function is purely logarithmic, this value corresponds 
to 21 — k powers of logarithms. In general, the space of possible transcendental functions is 
not fully characterized mathematically, particularly for functions of multiple dimensionless 
arguments. At the end of this subsection we give some examples of functions of definite 
transcendental weight, which appear in scattering amplitudes for massless particles, and 
which therefore might be considered candidates from which to build solutions for A. 

Because T, Ts p and A are associated with the 1/e single pole, their maximal transcen- 
dentality is T max — 21 — 1. For J\f — 4 super-Yang-Mills theory, in every known instance 
the terms arising in this way are purely of this maximal transcendentality: there are no 
terms of lower transcendentality. This property is relevant also for non-supersymmetric 
massless gauge theories, particularly at three loops. Indeed, in any massless gauge theory 
the purely-gluonic web diagrams that we need to consider at three loops are the same as 
those arising in M = 4 super- Yang- Mills theory. We conclude that at three loops A should 
have transcendentality r = 5 [58], while for I > 3 some relevant webs may depend on the 
matter content of the theory, so that A is only constrained to have a transcendentality at 
most equal to 21 — 1. 

It should be emphasized that some transcendentality could be attributed to constant 
prefactors. For example, the sum-over-dipoles formula (3.4) for r^p attains transcenden- 
tality r = 21 — 1 as the sum of r = 21 — 2 from the (constant) cusp anomalous dimension 7^ 
(associated with a 1/e 2 double pole in the amplitude) and r = 1 from the single logarithm. 

Because the functions A| m are defined up to possible numerical prefactors, which may 
carry transcendentality, terms of the form (5.7), (5.11) or (5.12) must obey 

h tot = h 1 + h 2 + h 3 <2l-l. (5.16) 

We note furthermore that constants of transcendentality r = 1, i.e. single factors of ir, 
do not arise in Feynman diagram calculations, except for imaginary parts associated with 
unitarity phases. We conclude that whenever the maximal transcendentality argument 
applies to T, the special case in which our functions A4 have r = 21 — 2 is not allowed. 

The sum of the powers of all the logarithms in the product must then be no more than 
^tot = 5 at three loops, or h tot = 7 at four loops, and so on. In the special cases considered 
above, at three loops, the constraint is: 3/i < 5, i.e. h < 1 in eq. (5.8), h < 5 in eq. (5.9), 
and 2h < 5, i.e. h < 2 in eq. (5.10). Clearly, at low orders, transcendentality imposes strict 
limitations on the admissible functional forms. We will take advantage of these limitations 
at three loops in Sec. 5.5. 

We close this subsection by providing some examples of possible transcendental func- 
tions that might enter A, beyond the purely logarithmic examples we have focused on so 
far. For functions of multiple dimensionless arguments, the space of possibilities is not pre- 
cisely characterized. Even for kinematical constants, the allowed structures are somewhat 
empirically based: the cusp anomalous dimension, for example, can be expressed through 
three loops [54] in terms of linear combinations of the Riemann zeta values C( n ) (having 
transcendentality n), multiplied by rational numbers; other transcendent als that might be 
present — such as In 2, which does appear in heavy-quark mass shifts — are not. 
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The cusp anomalous dimension governs the leading behavior of the twist-two anoma- 
lous dimensions for infinite Mellin moment N. At finite N, these anomalous dimensions can 
be expressed [54] in terms of the harmonic sums Sfi T (N) [68], where n T is a r-dimensional 
vector of integers. Harmonic sums are the Mellin transforms of harmonic polylogarithms 
Hm. T (x) [69], which are generalizations of the ordinary polylogarithms Li n (x). They are 
defined recursively by integration, 

R At (x) = f dx' f{a-x')^rh T _,{x'), (5.17) 

J 

where a = —1, or 1, and 

/(-l;x) = — |— , f(0; X ) = -, /(l ;a; ) = — !— . (5.18) 
1 + X x 1 — X 

Note that the transcendentality increases by one unit for each integration. All three values 
of a are needed to describe the twist-two anomalous dimensions. However, for the four- 
point scattering amplitude, which is a function of the single dimensionless ratio r defined 
in eq. (A.l), only a = 0, 1 seem to be required [70]. 

Scattering amplitudes depending on two dimensionless ratios can often be expressed in 
terms of harmonic polylogarithms as well, but where the parameter a becomes a function 
of the second dimensionless ratio [71]. In Ref. [72], a quantity appearing in a six-point 
scattering amplitude at two loops was recently expressed in terms of the closely-related 
Goncharov polylogarithms [73] in two variables, and at weight (trancendentality) four. 
Other recent works focusing more on general mathematical properties include Refs. [74,75]. 
In general, the space of possible functions becomes quite large already at weight five, and 
our examples below are meant to be illustrative rather than exhaustive. 

5.4 Collinear limits 

Equipped with the knowledge of how Bose symmetry and other requirements may be 
satisfied, let us return to the splitting amplitude constraint, namely the requirement that 
the difference between the two A terms in eq. (4.14) must conspire to depend only on the 
color and kinematic variables of the two partons that become collinear. 

We begin by analyzing the case of an amplitude with precisely four colored partons, 
possibly accompanied by other colorless particles (we postpone the generalization to an 
arbitrary number of partons to Sec. 6). The collinear constraint simplifies for n = 4 
because for three partons there are no contributions beyond the sum-over-dipoles formula, 
so that A n _i = A3 = [42]. 5 In eq. (4.14) we therefore have to consider A4 on its own, and 
require that when, say, p\ and P2 become collinear A4 does not involve the kinematic or 
color variables of other hard particles in the process. Because in this limit there remains no 
non-singular Lorentz-invariant kinematic variable upon which A4 can depend, it essentially 
means that A4 must become trivial in this limit, although it does not imply, of course, that 
A4 vanishes away from the limit. In the following we shall see how this can be realized. 

5 For n = 4 we should add a colorless particle carrying off momentum; otherwise the three-parton 
kinematics are ill-defined (for real momenta), and the limit pi ■ P2 — ► is not really a collinear limit but a 
forward or backward scattering limit. 
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To this end let us first carefully examine the limit under consideration. We work with 
strictly massless hard partons, pf = for all i. In a fixed- angle scattering amplitude we 
usually consider 2pi -pj = Q 2 (3i ■ (3j where Q 2 is taken large, keeping • fjj = 0(1) for any i 
and j. Now we relax the fixed-angle limit for the pair of hard partons p\ and p<i- Defining 
P = pi +P2 as in Sec. 4, we consider the limit 2pi -p^lQ 2 = P 2 /Q 2 — > 0. The other Lorentz 
invariants all remain large; in particular for any j ^ 1,2 we still have 2p\ ■ pj = Q 2 fii • (3 



j 



and 2p2 ■ Pj = Q 2 $2 • Pj where f3\ • f3j and /?2 ■ Pj are of O(l). In order to control the way 
in which the limit is approached, it is useful to define 



Pl = zP + k, p 2 = (l-z)P-k, (5.19) 

so that z measures the longitudinal momentum fraction carried by p\ in P, namely 

where we assume, for simplicity, that the "+" light-cone direction 6 is defined by p±, so that 
Pi = ipt i _ i 0j_). In eq. (5.20) both the numerator and denominator are of order Q, so z 
is of 0(1) and remains fixed in the limit P 2 /Q 2 ->• 0. In eq. (5.19) k is a small residual 
momentum, making it possible for P to be off the light-cone while p\ and P2 remain strictly 
light-like. Using the mass-shell conditions p\ = p\ = one easily finds 



k 2 = -z(l - z)P 2 , k-P= i(l-2z)P 2 , (5.21) 



so that the components of k are 



+ ,-^p T ,-V4l-z)P 2 ) . (5.22) 



Note that in the collinear limit k /Q scales as P 2 /Q 2 , while k±/Q scales as \J P 2 /Q 2 . 

We can now examine the behavior of the logarithms of the three cross ratios entering 
A4 111 in eq. (5.7), in the limit P 2 — > 0. Clearly, L1234 and L1423, which contain the vanishing 
invariant p\ -pi either in the numerator or in the denominator, will be singular in this limit. 
Similarly, it is easy to see that L1342 must vanish, because P1342 — > 1. More precisely, the 



6 One can then further specify the frame by choosing the "— " direction along the momentum of one of 
the other hard partons, say p$. 
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collinear behavior may be expressed using the parametrization of eq. (5.19), with the result 



-^1234 — hi 



Pi -P2P3- Pa 



yPl -P3P2-P4, 

~lnf P2 Pz-P* V k -P* + k -l± > oo (5.23) 

\2z(l- z)P -p 3 P -pa) zP-p 3 (1 - z) P ■ p A 



Lu23 = In 



Pi -PAP2-PZ 



Pi -P2P4-P3, 

-z) P- 

P 2 P4'P3 J (l-z)P-p 3 ZP-P4 



^( 2z(l-z)P-p 4 P-p 3 ^ _ k-p 3 | ky 4 (5.24) 



?3 PA ■ P2 \ _ 1 / k • p 3 fc • p 4 

P1-PAP3-P2J Z(l-Z) \P-p 3 P-PA 



V »i ■ »4 Da ■ »9 / z( 1 — Z ) V -r • p a P ■ Va I V / 

(5.25) 



where we expanded in the small momentum k. As expected, two of the cross-ratio loga- 
rithms diverge logarithmically with P 2 /Q 2 , with opposite signs, while the third cross-ratio 
logarithm vanishes linearly with ^/P 2 /Q 2 . We emphasize that this vanishing is indepen- 
dent of whether the momenta Pj are incoming or outgoing, except, of course, that the 
two collinear partons p\ and p2 must either be both incoming or both outgoing. Indeed, 
according to eq. (2.11), £1342 carries no phase when p\ and P2 are collinear: 



P1342 



Pi -P3PA- P2 



Pi -PAP3- P2 



g— i7r(Ai3+A42— A14— A32) ^ (5 26) 



since A13 = A32 and A42 = A14. 

Let us now examine a generic term with a kinematic dependence of the form (5.7) in 
this limit. Substituting eqs. (5.23) through (5.25) into eq. (5.7) we see that, if h 3 (the power 
of L1342) is greater than or equal to 1, then the result for A4" 1 in the collinear limit is zero. 
This vanishing is not affected by the powers of the other logarithms, because they diverge 
only logarithmically as P 2 /Q 2 — > 0, while L1342 vanishes as a power law in the same limit. 
In contrast, if h 3 = 0, and h\ or J12 is greater than zero, then the kinematic function A4" 1 
in eq. (5.7) diverges when p\ and P2 become collinear, due to the behavior of -L1234 and 
Li423 in eqs. (5.23) and (5.24). The first term in each of these equations introduces explicit 
dependence on the non-collinear parton momenta p 3 and p 4 into A n in eq. (4.14), which 
would violate collinear universality. We conclude that consistency with the limit where p\ 
and P2 become collinear requires h 3 > 1. 

Obviously we can consider, in a similar way, the limits where other pairs of partons 
become collinear, leading to the conclusion that all three logarithms, L1234, £1423 and L1324 
must appear raised to the first or higher power. Collinear limits thus constrain the powers 
of the logarithms by imposing 

hi>l, Vi. (5.27) 
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This result puts a lower bound on the transcendentality of A4 111 , 



namely 



hot = hi + h 2 + h 3 > 3 . 



(5.28) 



5.5 Three-loop analysis 

We have seen that corrections to the sum-over-dipoles formula involving four-parton cor- 
relations are severely constrained. We can now examine specific structures that may arise 
at a given loop order /, beginning with the first nontrivial possibility, / = 3. Because we 
consider webs that are attached to four hard partons, at three loops they can only attach 
once to each eikonal line, as in fig. 1, giving the color factor in eq. (5.2), where h abcd must 
be constructed out of the structure constants f ade . The only possibility is terms of the form 
jade jbce — same form we obtained in the previous section starting from the symmetry 
properties of the kinematic factors depending on L^^. In contrast, the symmetric tensor 
d ade cannot arise in three-loop webs. 

Taking into account the splitting amplitude constraint (5.27) on the one hand, and the 
maximal transcendentality constraint (5.16) on the other, there are just a few possibilities 
for the various powers h^. These are summarized in Table 1. 

The lowest allowed transcendentality for A| m is r = 3, corresponding to hi = h 2 = 
/13 = 1. This brings us to eq. (5.8), in which we would have to construct a completely 
antisymmetric tensor h abcd out of the structure constants f ade . Such a tensor, however, 
does not exist. Indeed, starting with the general expression (5.11), which is written in 
terms of the structure constants, and substituting hi = h 2 = h% = 1, we immediately 
see that the color structure factorizes, and vanishes by the Jacobi identity (5.13). The 
possibility hi = h 2 = /13 = 1 is thus excluded by Bose symmetry. 

Next, we may consider transcendentality r = 4. Ultimately, we exclude functions with 
this degree of transcendentality at three loops, because we are dealing with purely gluonic 
webs, which are the same as in N = 4 super Yang-Mills theory. We expect then that the 
anomalous dimension matrix will have a uniform degree of transcendentality r = 5, and 
there are no constants with r = 1 that might multiply functions with /i to t = 4 to achieve 
the desired result, as discussed in Sec. 5.3. However, it is instructive to note that symmetry 
alone does not rule out this possibility. Indeed, having excluded eq. (5.12), involving the 
symmetric tensor d ade , we may consider eq. (5.11) with hi + h 2 + h% = 4. Bose symmetry 
and the splitting amplitude constraint in eq. (5.27) leave just two potential structures, one 
with hi = 2 and h 2 = h$ = 1, and a second one with h\ = h 2 = 1 and /13 = 2 (h\ = /13 = 1 
and h 2 = 2 yields the latter structure again). The former vanishes identically, while the 
latter could provide a viable candidate, 



Ai m) (p^) = T?T b 2 T^ 




(5.29) 
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hi 


h 2 


h 3 


htot 


comment 


1 

o 
Z 

1 

1 

3 
2 
1 


1 
1 
1 
2 
1 
1 
1 


1 
l 

2 
2 
1 
2 
3 


3 
4 
4 
5 
5 
5 
5 


vanishes identically by Jacobi identity (5.13) 

kinematic factor vanishes identically 

allowed by symmetry, excluded by transcendentality 

viable possibility, eq. (5.30) "*\ 

viable possibility, eq. (5.31) 1 . . . 

. , , (K oon > a11 coincide using eq. (5.13) 
viable possibility, eq. (5.33) g 

viable possibility, eq. (5.34) J 



Table 1: Different possible assignments of the powers hi in eq. (5.11) at three loops. We only 
consider hi > 1 because of the splitting amplitude constraint (5.27) and h to t < 5 because of the 
bound on transcendentality, eq. (5.16). We also omit the combinations that can be obtained by 
interchanging the values of hi and hs; this interchange yields the same function, up to a possible 
overall minus sign. 



We rule out eq. (5.29) based only on its degree of transcendentality. 

We consider next the highest attainable transcendentality at three loops, r = 5. Equa- 
tion (5.11) yields four different structures, summarized in Table 1. The first structure we 
consider has h± = 1 and h 2 = h$ = 2. It is given by 



A 



(122) 



(Pijkl) 



i fcae fdbe t it t \l , fbae fcae t it t \ 

+ J J -^1423 1-^1234 -^1342 J +J J -^1342 (,-^1423 -^1234 J 

The second structure has h\ = 3 and h 2 = /13 = 1, yielding 



rjia rj-16 rj-iC rj-iG? 



/ade fcbe j it r \^ 

/ -^1234 (,-^1423 -^1342 ) 



pbae fcde - 



(5.30) 



A 



(311) 



(Pijkl) 



rjia r£i6 rj-iC rj-iC? 



jade J c b e ^]^ 



1234) -^1423-^1342 



1 rcae fdbe / t \6 t t . fOaefcaeij \i r r 

+ J J 1,-^1423 ) -^1234 -^1342 + / J (,-^1342 J -^1423-^1234 



6ae j^cde / 



(5.31) 



We now observe that the two functions (5.30) and (5.31) are, in fact, one and the 
same. To show this, we form their difference, and use relation (5.6) to substitute L1234 = 
—-^1423 — £1342 ■ We obtain 



A 



(122) 



A 



(311) 



TCLrrib rriCrrid T T T 

1 1 2 1 3 1 4 -^1234 -^1423 -^1342 



X 



+ / uu,c /" ic [ L 142 3 L1234 — L 



fade rcbe ( j j j 2 \ . rcae fdbe I j j j 2 \ 

J J 1-^1423 -^1342 — -^1234 ) + / / 1-^1234 -^1342 ~ -^1423 J 

jbae jcde ^ j 

Tarribrricrrid t t t fade fcbe i fcae fdbe , fbae fcde 

1 1 2- L 3- L 4 -^1234-^1423-^1342 / / +J J +J J 

x (-^1342 + -^1342 -^1423 + -^1423) = > 



2 ' 
1342 



(5.32) 



vanishing by the Jacobi identity (5.13). 



- 24 - 



The last two structures in Table 1 are given by 



Al 212 W0 = T?T b 2 T^ 



arr\brr\crr\d 



fade i-cbe j 2 
J J -^1234 

cae rdbe t 2 



. (-^1423 -^1342 + -^1423 ^1342^ 
+ / Cae / a ° e -^1423 (-^1234 -^1342 + ^1234 -^1342 
+ f baC f Cde L\ M2 (-^1423 -^1234 + -^1423 L 12Zij 



(5.33) 



and 



^ 13 \pijkl) — T^T^TgT^ f ade f Cbe Li234 (-^1423 -^1342 + -^1423-^1342^) 

+ f Cae f dbe Lu23 (-^1234 -^1342 + -^1234 ^1342^ 



(5.34) 



_l_ jbae jcdej^ 



1342 I -^1423 -^1234 



+ L 



1423 -^1234 



One easily verifies that they are both proportional to A^ 122 -* = A^ U \ Consider first 
eq. (5.33). In each line we can factor out the logarithms and use eq. (5.6) to obtain a 
monomial. For example, the first line may be written as: 



2/22 \ 2 / \ 

-^1234 ( -^1423 -^1342 + -^1423 -^1342 J = -^1234 -^1423 -^1342 (-^1423 + -^1342) 



"-^1234 -^1423 -^1342 



(5.35) 



where we recognise that this function coincides with eq. (5.31) above. Consider next 
eq. (5.34), where, for example, the first line yields 

-^1234 ^1423 -^1342 + -^1423 -^1342 ) = ^1234 -^1423 -^1342 ^1342 + ^1423 ) 

= Ll234 -^1423 -^1342 ( (-^1342 + -^1423 ) 2 ~ 2L1342 ^1423 
= L1234 -^1423 -^1342 ^1234 ~~ 2L1342 ^1423^) , 



(5.36) 



which is a linear combination of eqs. (5.30) and (5.31), rather than a new structure. 



We conclude that there is precisely one function, A 



(122) 



A? n \ that can be con- 



structed out of arbitrary powers of logarithms and is consistent with all available constraints 
at three loops. We emphasize that this function is built with color and kinematic factors 
that one expects to find in the actual web diagram computations, and it is quite possible 
that it indeed appears. Because this structure saturates the transcendentality bound, its 
coefficient is necessarily a rational number. 

Note that color conservation has not been imposed here, but it is implicitly assumed 
that for a four-parton amplitude 



+ T 2 + T3 + T4 = . 



(5.37) 



Importantly, upon using this relation, the structure (5.30) (or, equivalently, (5.31)) remains 
non-trivial. 
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5.6 Three-loop functions involving polylogarithms 

Additional functions can be constructed upon removing the requirement that the kinematic 
dependence be of the form (5.7), where only powers of logarithms are allowed. Three key 
features of the function In p were essential in the examples above: it vanishes like a power 
at p = 1, it has a definite symmetry under p — > 1/p, and it only diverges logarithmically 
as p — > and p — > oo. These properties can be mimicked by a larger class of functions. In 
particular, allowing dilogarithms one can easily construct a function of transcendentality 
r = 4, which is consistent with Bose symmetry and collinear constraints. It is given by 

Af 1,Li2 W) = T?T^Tf 3 Tf (5.38) 

f ade f Cb£ (Li2(l " P1234) " Li 2 (l - 1/P1234)) lnpi423 In P1342 
+ f Cae f dbe (Li 2 (l - p 14 23) " Li 2 (l - 1/P1423)) ln/>1234 In P1342 
+/ &a 7 Cde (Li 2 (l-pl342)-Li 2 (l-l/pi342)) lnpi423 In P1234 • 

The key point here is that the function Li 2 (l — pi 2 34) — Li 2 (l — 1 //O1234) is odd under 
P1234 — > I/P1234, which allows it to be paired with the antisymmetric structure constants 
f ade . It is also easy to verify that the collinear constraints are satisfied. 

We note that it is also possible to construct a potentially relevant function containing 
logarithms with a more complicated kinematic dependence. Indeed, the structure 

Ai 211 ' mod) (p ijjH ) = T?T fe 2 T^ (5.39) 
f ade f cbe In p 1234 ln ( (13^)2) ln ^423 lnp 1342 



+f cae f dbe lnpu23 tn ( ^ j 
A 1 -Pl423j 



+f bae f cde lnp 1U2 ln( P1U2 ) lnp 142 3 lnpi 23 4 
V(l -P1342) 2 / 

fulfills the symmetry requirements discussed above, because In (pi234/(l — P1234) 2 ) is even 
under P1234 — > l/pi234- Thanks to the extra power of the cross-ratio logarithm, eq. (5.39) 
also vanishes in all collinear limits, as required. Logarithms with argument 1 — p^^i cannot 
be directly rejected on the basis of the fact that they induce unphysical singularities, 
because ptjki — > 1 corresponds to a physical collinear limit 7 . We conclude that eqs. (5.38) 
and (5.39) would be viable based on symmetry and collinear requirements alone. However, 
we can exclude them on the basis of transcendentality: as discussed in Sec. 5.3, a function 
with h tot = 4 cannot arise at three loops, because it cannot be upgraded to maximal 
transcendentality r = 5 by constant prefactors. 



7 Note that the analogous structure containing In (pi234/(l + P1234) 2 ) can be excluded because the limit 
pijki — > — 1 should not be singular. Indeed, by construction the variables always contain an even 

number of negative momentum invariants, so their real part is always positive (although unitarity phases 
may add up and bring their logarithm to the second Riemann sheet). 
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At transcendentality r = 5, there are at least two further viable structures that involve 
poly logarithms, in which second and third powers of logarithms are replaced, respectively, 
by appropriate combinations Li2 and Li3. The first structure can be obtained starting from 
cq. (5.30), and using the same combination of dilogarithms that was employed in eq. (5.38). 
One finds 

Ai 122>Li2 W,) = T1T2T3T4 (5.40) 

f ade f cbe lnpi234 _ ^ _ ^ _ y ^(1 - p U23 ) - Li 2 (l - l/pi 4 2 3 )) 

+ f cae f dbe _ ^ _ ^ _ y^fj (^(1 _ _ 1^(1 _ l/p 1342 )) 

+ / fea 7 Cde lnpi342 (Li 2 (l " P1234) " Li 2 (l - 1/P1234)) (u 2 (l ~ Pl42 3 ) ~ Li 2 (l - l/pi42 3 )) 

Here it was essential to replace both In 2 terms in order to keep the symmetry properties in 
place. Starting instead from eq. (5.31), there is one possible polylogarithmic replacement, 
which, however, requires introducing trilogarithms, because using Li 2 times a logarithm 
would turn the odd function into an even one, which is excluded. One may write instead 

Al 311 ' Li3 W<) = T1T2T3T4 (5.41) 

X f ade f Cbe (Li 3 (l - p 123 4) " Li 3 (l " 1/P1234)) ^1423 W 
+ / Cae / d6e (Li 3 (l " P1423) " Li 3 (l " 1/P1423)) ^1234 ^1342 
+ / fea 7 Cde (Li 3 (l " P1342) " Li 3 (l - 1/P1342)) ^1423 ^1234 

Neither eq. (5.40) nor eq. (5.41) can be excluded at present, as they satisfy all available 
constraints. We can, however, exclude similar constructions with higher-order polyloga- 
rithms. For example, L14 has transcendentality r = 4, so it could be accompanied by at 
most one logarithm; this product would not satisfy all collinear constraints. We do not 
claim to be exhaustive in our investigation of polylogarithmic functions; additional pos- 
sibilities may arise upon allowing arguments of the polylogarithms that have a different 
functional dependence on the cross ratios. 

5.7 Four-loop analysis 

Let us briefly turn our attention to contributions that may arise beyond three loops. At the 
four-loop level several new possibilities open up. First, there are potential quartic Casimirs 
in jk- Corresponding corrections to the soft anomalous dimension would satisfy inhomo- 
geneous differential equations, eq. (5.5) of Ref. [42]. Beyond that, new types of corrections 
may appear even if jk admits Casimir scaling. First, considering the logarithmic expres- 
sions of eq. (5.7), purely gluonic webs might give rise to functions of transcendentality up 
to /itot = h± + /t 2 + hs = 7. At this level, there are four potential functions: (a) h± = 5 
and /i 2 = /13 = 1; (b) hi = 4, hi = 2, /13 = 1; (c) h\ = 1 and hi = /13 = 3; (d) hi = 3 
and hi = h^ = 2. Of course, as in the three-loop case, also polylogarithmic structures may 
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appear, and functions with /i to t < 5 might be present (of the type already discussed at 
three loops), multiplied by transcendental constants with r > 2. 

It is interesting to focus in particular on color structures that are related to quartic 
Casimir operators, which can appear at four loops not only in -jk but also in four-parton 
correlations. Indeed, a structure allowed by Bose symmetry and collinear constraints is 
given by eq. (5.8), where the group theory factor h abcd is generated by a pure-gluon box 
diagram attached to four different hard partons, giving rise to a trace of four adjoint 
matrices. It is of the form 

A± 4 ' A (p ljM ) = Tr [F a F b F c F d ] T?T 6 2 T^ [lnp 1234 In p U23 lnp 13 4 2 ]\ (5.42) 

where F a are the SU(iV c ) generators in the adjoint representation, (F a )t, c = —if a bc - This 
expression may be relevant a priori for both odd and even h, projecting respectively on the 
totally antisymmetric or symmetric parts of Tr [F a F b F c F d ] . As noted above, however, 
a totally antisymmetric tensor cannot be constructed with four adjoint indices, so we are 
left with the completely symmetric possibility, which indeed corresponds to the quartic 
Casimir operator. The transcendentality constraint comes into play here: the only even 
integer h that can give transcendentality r < 7 is h = 2. For N = 4 super- Yang-Mills 
theory, eq. (5.42) with h = 2 can be excluded at four loops, because there is no numerical 
constant with r = 1 that could bring the transcendentality of eq. (5.42) from 6 up to 7. 

On the other hand, in theories with a lower number of super symmetries, and at four 
loops, there are potentially both pure-glue and matter-loop contributions of lower tran- 
scendentality, because only the specific loop content of N = 4 super- Yang-Mills theory 
is expected to be purely of maximal transcendentality. Thus eq. (5.42) may be allowed 
for h = 2 for generic adjoint-loop contributions (for example a gluon box in QCD), and 
analogously 



A^' F (p ijkl ) = Tr \t a t b t c t d \ T?T 6 2 T^ [In p 1234 lnp U23 hipi 3 4 2 J , (5.43) 

for loops of matter in the fundamental representation (with generators t a ), e.g. from quark 
box diagrams. As before, the other power allowed by transcendentality, h = 1, is excluded 
by symmetry, because there is no projection of Tr [t a t b t c t rf ] that is totally antisymmetric 
under permutations. 

While eq. (5.42) is excluded by transcendentality for N = 4 super- Yang-Mills theory, 
another construction involving the quartic Casimir is allowed: eq. (5.15) for h\ = 2, h 2 = 
/13 = 1 can be used, after multiplying it by the transcendentality r = 3 constant C(3)- 
Finally, as already mentioned, there are a number of other purely logarithmic structures 
with partial symmetry in each term, as represented by eqs. (5.11) and (5.12), that may 
appear at four loops. 

6. Generalization to n-leg amplitudes 

The above analysis focused on the case of four partons, because this is the first case where 
cross ratios can be formed, and thus A may appear. However, rescaling-invariant ratios 



h 
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involving more than four momenta can always be split into products of cross ratios involving 
just four momenta. Therefore it is straightforward to generalize the results we obtained to 
any n-parton process at three loops. Indeed, contributions to A n are simply constructed 
as a sum over all possible sets of four partons, 

A n = MPijkl), (6.1) 

just as the sum-over-dipoles formula (2.7) is written as a sum over all possible pairs of legs. 
The indices in the sum in eq. (6.1) are of course all unequal. Assuming a purely logarithmic 
structure, at three loops the function A4 in eq. (6.1) is given by A 4 in eq. (5.30) (or, 
equivalently, A^ 311 '* in eq. (5.31)). Of course the overall prefactor to A4 could still be zero; 
its value remains to be determined by an explicit computation. The total number of terms 
in the sum increases rapidly with the number of legs: for n partons, there are (4) different 
terms. 

Now we wish to show that this generalization is a consistent one. To do so, we shall 
verify that the difference between A n and A n _i in eq. (4.14), for the splitting amplitude 
anomalous dimension, does not introduce any dependence on the kinematics or color of 
any partons other than the collinear pair. The verification is non-trivial for n > 5 because 
A n _i is no longer zero. 

Consider the general n-leg amplitude, in which the two legs p\ and P2 become collinear. 
The terms entering eq. (4.14) include: 

• A sum over ( n ^ 2 ) different terms in A n that do not involve any of the legs that 
become collinear. They depend on the cross ratios Pijki where none of the indices is 1 
or 2. However, exactly the same terms appear in A n _i, so they cancel in eq. (4.14). 

• A sum over { n ' 2 2 ) different terms in A n depending on the variables p\2ij (and per- 
mutations), where i, j 7^ 1, 2. These variables involve the two legs that become 
collinear. According to eq. (6.1), each of these terms is A4, namely it is given by a 
sum of terms that admit the constraint (5.27). Therefore each of them is guaranteed 
to vanish in the collinear limit, and we can discard them from eq. (4.14). The same 
argument applies to any A4 that is consistent at the four-parton level, such as the 
polylogarithmic constructions (5.40) and (5.41). 

• Finally, A n brings a sum over 2 x ( n ^ 2 ) terms involving just one leg among the two 
that become collinear. These terms depend on pijki or p2jkh where j, k, I 7^ 1, 2. In 
contrast, A n _i brings just one set of such terms, because the (12) leg, P = p\ + P2, 
is now counted once. Recalling, however, that this leg carries the color charge 

T a = T" + T2 , (6.2) 

it becomes clear that any term of this sort having a color factor of the form (5.2) 
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would cancel out in the difference. Indeed 

A n (PljkhP2jkl) - A n _i (p(i2)jkl) 



EJ^abcd rpbrpc rj-i(i rjlO 
j k I ^ '■ 



tA^(p ljkl )+T^A^(p2 jkl ) - T a A^(p {12)jkl ) (63) 



j,k,l 




= 0. 



To show that this combination vanishes we used eq. (6.2) and the fact that the 
kinematic factor A km in all three terms is identical because of rescaling invariance, 
that is, it depends only on the directions of the partons, which coincide in the collinear 
limit. 

We conclude that eq. (6.1) is consistent with the limit as any two of the n legs become 
collinear. 

A similar analysis also suggests that eq. (6.1) is consistent with the triple collinear 
limit in which pi, p2 and p% all become parallel. We briefly sketch the analysis. We assume 
that there is a universal factorization in this limit, in which the analog of Sp again only 
depends on the triple-collinear variables: P 2 = {pi + P2 + P3) 2 , which vanishes in the limit; 
2pi • P2/P 2 and 2p 2 -p^/P 2 ] and the two independent longitudinal momentum fractions for 
the pi, namely z\ and Z2 (and Z3 = 1 — z\ — 22) — see e.g. Ref. [76] for a discussion at one 
loop. In the triple-collinear limit there are the following types of contributions: 

• (™4 3 ) terms in A n that do not involve any of the collinear legs. They cancel in the 
analog of eq. (4.14) between A n and A n _2, exactly as in the double-collinear case. 

• 3 x ( n g 3 ) terms containing cross ratios of the form pijki, or similar terms with 1 
replaced by 2 or 3. These contributions cancel exactly as in eq. (6.3), except that 
there are three terms and the color conservation equation is T a = Tf + Tf + T|. 

• 3x( n ~ 3 ) terms containing cross ratios of the form p\2kh or similar terms with {1, 2, 3} 
permuted. These terms cancel for the same reason as the pyiij terms in the double- 
collinear analysis, namely one of the logarithms is guaranteed to vanish. 

• [n— 3) terms containing cross ratios of the form p\2zi- This case is non-trivial, because 
no logarithm vanishes (no cross ratio goes to 1). However, it is easy to verify that in 
the limit, each of the cross ratios that appears depends only on the triple-collinear 
kinematic variables, and in a way that is independent of pi. Therefore the color 
identity 2 3 ^1 = ~~ T a can ^ e use< ^ to express the color, as well as the kinematic 
dependence, of the limit of A n solely in terms of the collinear variables, as required 
by universality. 

Thus all four contributions are consistent with a universal triple-collinear limit. However, 
because the last type of contribution is non-vanishing in the limit, in contrast to the double- 
collinear case, the existence of a non-trivial A n would imply a new type of contribution to 
the 1/e pole in the triple-collinear splitting function, beyond that implied by the sum-over- 
dipoles formula. 
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We conclude that eq. (6.1) provides a straightforward and consistent generalization 
of the structures found in the four-parton case to n partons. At the three-loop level, if 
four-parton correlations arise, they contribute to the anomalous dimension matrix through 
a sum over color 'quadruples' of the form (6.1). At higher loops, of course, structures 
that directly correlate the colors and momenta of more than four partons may also arise. 

7. Conclusions 

Building upon the factorization properties of massless scattering amplitudes in the soft and 
collinear limits, recent work [42,43] determined the principal structure of soft singularities in 
multi-leg amplitudes. It is now established that the cusp anomalous dimension 7^ controls 
all pairwise interactions amongst the hard partons, to all loops, and for general N c . The 
corresponding contribution to the soft anomalous dimension takes the elegant form of a 
sum over color dipoles, directly correlating color and kinematic degrees of freedom. This 
recent work also led to strong constraints on any additional singularities that may arise, 
thus opening a range of interesting questions. 

In the present paper we studied multiple constraints on the form of potential soft 
singularities that couple directly four hard partons, which may arise at three loops and 
beyond. We focused on potential corrections to the sum-over-dipoles formula that do not 
require the presence of higher Casimir contributions to the cusp anomalous dimension jk- 
The basic property of these functions is that they satisfy the homogeneous set of differential 
equations (2.9), and therefore they can be written in terms of conformally invariant cross 
ratios [42]. 

Our main conclusion is that indeed, potential structures of this kind may arise starting 
at three loops. Their functional dependence on both color and kinematic variables is, 
however, severely constrained by 

• Bose symmetry; 

• Sudakov factorization and momentum-rescaling symmetry, dictating that corrections 
must be functions of conformally invariant cross ratios; 

• collinear limits, in which the (expected) universal properties of the splitting amplitude 
force corrections to vanish (for n = 4 partons) or be smooth (for n > 4 partons) in 
these limits; 

• transcendentality, a bound on which is expected to be saturated at three loops, based 
on the properties of N = 4 super- Yang-Mills theory. 

In the three-loop case, assuming purely logarithmic dependence on the cross ratios, these 
constraints combine to exclude all but one specific structure. The three-loop result for A n 

('122) 

can therefore be written in terms of the expression A 4 in eq. (5.30), up to an overall 
numerical coefficient. Because this structure has the maximal possible transcendentality, 
r = 5, its coefficient is a rational number. For all we know now, however, this coefficient 
may vanish. It remains for future work to decide whether this contribution is present or 
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not. Considering also poly logarithmic functions of conformally invariant cross ratios in 
A4, we find that at three loops at least two additional acceptable functional forms arise, 
eqs. (5.40) and (5.41). 

The range of admissible functions at four loops is even larger. A particularly inter- 
esting feature at this order is the possible appearance of contributions proportional to 
quartic Casimir operators, not only in the cusp anomalous dimension, but in four-parton 
correlations as well. 

Explicit computations at three and four loops will probably be necessary to take the 
next steps toward a complete understanding of soft singularities in massless gauge theories. 
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A. The four-parton amplitude with no momentum recoil 

Here we investigate briefly the simplifications of the potential forms for A4 at three loops 
that result from using momentum-conservation relations special to the four-parton ampli- 
tude, pi + P2 + P3 + Pa = 0. Thus we exclude here the presence in the amplitude of other 
colorless particles that might carry recoil momentum. As before, all momenta are light-like, 
pf = 0, so the momentum invariants are now related by 

P1-P2+P2-P3+P1-P3 = , 

as well as P3 ■ P4 = Pi ■ P2, etc. Using this relation, all three cross ratios entering A4 can 
be expressed in terms of a single dimensionless ratio, 



P2 -P3 



Pi -P3 



(A.l) 



r = 



Pi -P2 



Pi -P2 



Substituting into eq. (2.11), we have 




(A.2) 
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where we used the variables A^, defined below eq. (2.3), to keep track of the unitarity 
phases. 
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Figure 3: The logarithmic functions A4 with transcendentality r = 5 introduced in eqs. (5.30) 
and (5.31), respectively, in the case of a four-parton amplitude with no recoil. In each case we 
separately display the real (dot-dash) and imaginary (solid line) parts as a function of the ratio r, 
defined in eq. (A.l), multiplying the color coefficient f ade f cbe after the Jacobi identity has been 
taken into account, as in eq. (A. 3). In each case symbols represent eq. (5.30) while the lines stand 
for eq. (5.31). The plot demonstrates that the two are identical once the Jacobi identity is taken 
into account. 

Let us now examine the behavior of the three expressions for A4 that we found admissi- 
ble, as functions of r. As mentioned in Sec. 5.4, in the four-parton case with no momentum 
recoil, pi ■ pj is not a collinear limit, but a forward or backward scattering limit. There are 
three channels to consider: 

a) pi and P2 incoming, p% and P4 outgoing. The physical region is — 1 < r < 0. 

b) p\ and p% incoming, P2 and P4 outgoing. The physical region is < r < 00. 

c) pi and P4 incoming, P2 and P3 outgoing. The physical region is —00 < r < — 1. 

The two endpoints of each physical interval are the forward and backward scattering limits. 
Using eq. (A. 2), we can read off the phases associated with each of the logarithms of the 
cross ratios in these three physical regions. The results are summarized in Table 2. 

We have seen that, due to the Jacobi identity, the three terms in A4 corresponding to 
antisymmetrization of any pair of color indices are related to each other, leaving just two 
independent terms. For example, in eq. (5.30) we can substitute f bae f cde = _ja.de jebe _ 
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a) -1 < r < 



b) < r < oo 



c) — oo < r < — 1 



£1234/2 
£1423/2 
£1342/2 



In |(1 + r)/r 



— In |1 + r\ — in 
In Irl + in 



— In |1 + r| + in 
In |r| 

In |(1 + r)/r| — i7r 



— In 1 1 + r\ 

In |r| — in 

In |(1 + r)/r| + in 



Table 2: Analytic continuation of the three cross ratios into the three physical regions in a 2 — > 2 
scattering amplitude. 



The resulting dependence on r is shown in Figure 3, which displays the coefficient of 
jadejcbe^ a ft er us i n g the Jacobi identity, as in eq. (A. 3). The plot shows that the result is 
the same when starting with either eq. (5.30) or eq. (5.31), as proven in (5.32). In a similar 
way one can use the Jacobi identity for the other admissible functions given in eqs. (5.40) 
and (5.41). Of course, each of them yields a different function of r. 
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